Statistics of Lyapunov exponent in one-dimensional layered systems 
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Localization of acoustic waves in a one dimensional water duct containing many randomly dis- 
tributed air filled blocks is studied. Both the Lyapunov exponent and its variance are computed. 
Their statistical properties are also explored extensively. The results reveal that in this system the 
single parameter scaling is generally inadequate no matter whether the frequency we consider is 
located in a pass band or in a band gap. This contradicts the earlier observations in an optical case. 
We compare the results with two optical cases and give a possible explanation of the origin of the 
different behaviors. 
PACS numbers: 42.25.Dd, 89.75.Da 



I. INTRODUCTION 

The fact that the electronic localization in disordered 
systems Q is of wave nature has led to suggestion that 
classical waves could be similarly localized in random sys- 
tems. The effort in searching for localization of classi- 
cal waves such as acoustic and electro-magnetic waves is 
tremendous. It has drawn intensive attentions from both 
theorists [|]-D an d experimentalists |S JTo[ ] . Snice the pi- 
oneering work of Anderson et.al. the concepts of 
universality and scaling have become important jl2] [tG| . 
These ideas stem from the insensitivity of the macro- 
scopic laws to microscopic details; that is, systems or 
models which differ from each other on a microscopic 
level can show identical macroscopic behavior. Accord- 
ing to the hypothesis of single-parameter scaling (SPS) 
pd| , p^|p6| , if the localization behaviors of a ID disordered 
system obey SPS, the Lyapunov exponent or inverse lo- 
calization that characterizes the degree of localization 
will be proportional to its variance. 

In a recent study of acoustic (AC) waves propagation in 
a one dimensional randomly layered system |l7|] we have 
found that the statistics relation between Lyapunov ex- 
ponent (LE) 7 and its variance var(7) do not follow the 
preditions of single parameter scaling (SPS). However, in 
an earlier study on ID localization behaviors of electro- 
magnetic (EM) waves Jll|, the author claimed that the 
non-universal bahaviors of LE will disappear and SPS 
will be restorted while the randomness of the system ex- 
ceeds a critical value. Since for ID propagation the AC 
and EM waves are in fact mathematically equivalent, i.e., 
there exists one to one correspondence between these two 
kinds of waves, it looks quite impossible that they can 
have different localization behaviors. In order to under- 
stand where do the main differences between these two 
kinds of models come from, in this paper we study EM 
and AC systems simultaneously. Two EM models and 
one AC model are studied in this paper. We find that 
though the statements made by Deych et.al. [|15[ are cor- 
rect in their chosen case, however, the applicability of 
SPS is more or less based on the fact that the impedance 



contrast between the constituents of the wave media is 
closed to 1. Without this restriction then even in the EM 
systems SPS will not be restored in the high randomness 
limit. 

This paper is organized as follows. In the next sec- 
tion we explain the correspondence between EM and AC 
waves and define the three models employed in this pa- 
per. In Sec. Ill we we first review the previous results 
and then discuss the numerical results of the three cho- 
sen models. A possible explanation of the origin of the 
novel properties of the AC localization are also proposed. 
Concluding remarks are given in section IV. 



II. THEORY AND MODELS 

To begin with, we explain the one to one correspon- 
dence of the ID propagation between AC and EM waves. 
For simplicity while without destroying the generality we 
restrict our discussion to the monochromatic waves with 
time dependence e~ luJt . We also assume that waves are 
nomally incident on the left boundary of the media and 
propagate along the x-axis. The ID propagation of AC 
waves under these assumptions is govened by 
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where p = p(x), p = p(x) and c = c{x) represent the mass 
density, the pressure, and the phase velocity of the wave 
in the media, respectively. In the special case of layered 
media considered in p and c are all constants in a 
single layer. Across an interface that separate two layers 
either p or c jump to a diferent value but pressure p and 
media vibration velocity 



1 dp 

iujp dx 



(2) 



must be continuously connected. The continuity condi- 
tions for every interfaces together with the wave equation 
(fy) itself determine the dynamics of the whole system. 
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Similarly, the equation governing ID propagation of 
monochromatic EM waves can be deduced from Maxwell 
equation and is written as 



d 
dx 



1 dE 

p dx 



u 2 eE 



m 2 E 

jJLC 2 
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In a single layer the impedance of AC waves is given 
by pc and the impedance of EM waves is given by ec or 
pc. Thus we see that in Model 3 the impedance contrast 
is about 3365 and impedance contrast in Model 1 and 2 
are about 1.414 and 4.472 respectively. 
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where E — E y and H — H z are the electric and magnetic 
fields , e and fi stand for the permittivity and permeabil- 
ity, Co is the speed of light in vacuum and c = c^j JTp is 
the speed of light in the media. Comparing Eq. (0) and 
(^) with Eq. (|l|) one can easily recognize the equivalence 
between AC and EM models via the substitutions 



1/pc 2 , p^p 



(5) 



E — > p, H — > —u, 
in Eq. (g) or 

H —> p, E — it, e -> p, fi-> l/pc 2 (6) 
in Eq. (g). 

In spite of these similarities, the two models discussed 
in 1 17 1 and Jll| indeed have some different features. First, 
in the optical case that described in Eq. (j^), one usually 
assumes p = I. However, in our AC model the corre- 
sponding quantity is p. The mass density ratio btween 
water and air is about 775, a very large value. Second, 
the phase velocity ratio in our model is 4.455, whereas in 
Jl5| the ratio is 1.0945. Third, the thickness ratio in our 
model is 9999/1, much larger than 1/1 that considered 
in Jl5| ] . Fourth, we randomize the thickness of the water 
layers (medium with high phase velocity) and keep the 
thickness of air layers (medium with low phase velocity) 
constant. In Mj they randomize the thickness of the 
layers with low phase velocity. 

To clarify where do the major differeces come from, 
we define three models and study them numerically. The 
system for each model is a composite made of two kinds of 
material A and B with corresponding thickness aj and bj 
in the jth A/B layer. For simplicity hereafter we assume 
dj = a, i.e., all A type layers have the same size. Any 
quantity Q in A and B type layers are denoted as Q a and 
Qb respectively. 

• Model 1 is an optical model with Eb/e a — 2 and 
Hb/n<a = 1- The thickness ratio is (b)/a = 1, where 
be (6)[1 — A, 1 + A] and A G [0, 1]. 

• Model 2 is also an optical model with e&/e a = 20 
and pb/pa = 1- The thickness ratio is (b)/a = 1, 
where b G (b)[l - A,l + A] and A G [0, 1]. 

• Model 3 is our previously considered acoustic model 
with Cb/c a — 4.455 and Pb/p a = 755.2. The 
thickness ratio is given by (b)/a = 9999, where 
be (£>)[! — A, 1 + A] and A G [0, 1]. 



III. NUMERICAL RESULTS 

Before the discussion of the numerical results we first 
summarize the relevant results of our previous study p7[ . 
There both LE and its variance as functions of frequency 
were studied. At low disorders, the variance of LE in- 
side the gaps is small. Contrast to the optical case [ jUsl , 
there are no double maxima inside the gap. With in- 
creasing disorder, double peaks appears inside the allow 
bands. When exceeding a certain critical value, how- 
ever, the double peaks emerge. The higher frequency, 
the lower is the critical value. The increasing disorder re- 
duces the band gap effect and smears LE. We also ploted 
LE- variance relations. However, with increasing disor- 
der, we did not observe linear dependence between LE 
and its variance, as expected from the single parameter 
scaling theory. 

Now we turn to the discussion of the three chosen mod- 
els. When randomness A = 0, the layered systems be- 
come periodic. Eigenfunctions of wave equation in a pe- 
riodic environment are Bloch waves. As is well known, 
band structure appears in this situation. The under- 
standing of band structure is very important and help- 
ful in the following discussions. The dispersion relation 
of the Bloch waves in the underlying periodic system is 
given by 

cos Kd = cos k a a cos kbb — cosh 2?y sin k a a sin kbb (7) 

with d = a + b representing the thickness of the space 
period and K the Bloch wave number. Here for model 1 
and 2 function cosh 2rj is defined as 



cosh 2r? = - ( . I— 
2 V V e, 



and for model 3 it is given by 




(8) 



cosh2?7 = i (ah + ^ J , g = p a j ' pb, h = c a /cb- (9) 

In the frequency ranges that Kd are real, the waves are 
freely propagating in the media and by definition the fre- 
quency ranges correspond to the pass bands. Beyond the 
pass bands Kd are not purely real numbers, solutions of 
wave equation that satisfy appropriate boundary condi- 
tions in both the left and right infinity do not exist, and 
thus the frequency ranges are refer to as the band gaps. 
If the media is semi-infinite and the periodicity is ended 
by a boundary, say, the left boundary, then the waves 
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are localized in the vicinity of the boundary. The pene- 
tration depth is equal to l/|Im(AT)|. Fig. l(al)-(cl) plot 
the band structures of model 1-3. Solid curves represent 
Ke(Kd) and cover the pass bands. Broken lines cover the 
band gaps and represent the inverse penetration depths 
lm(Kd). Model 1 has very wide pass bands and very 
narrow band gaps. Band gaps are wider than pass bands 
in model 2. Model 3 has very wide band gaps and very 
narrow pass bands. An important feature of model 3 is 
that in the band gaps the the penetration depth is very 
small. Even in the first gap (which has the longest pene- 
tration depth) the penetration depth is smaller than one 
period. 

To see how randomness influences the transmission 
properties we select the frequency range that around the 
second gap and plot the transmission curves on (a2)-(c2). 
In this calculation 1600 layers (800 periods) are used in 
model 1 and 200 layers (100 periods) are used in model 2 
and 3. In the gaps the transmission rate is almost zero. 
When randomness is increased we observe that the trans- 
mission rate in the pass bands have been reduced much 
in the model 2 and model 3. On the other hand, the band 
gaps in model 2 and 3 seem to be more robust than in 
model 1. 

To further explore the influence of randomness we 
study both LE and its variance. The results are shown in 
Fig. g. (A1)-(A5) are results of model 1 and (B1)-(B5) 
are results of model 2. LE is denoted as 7 and its vari- 
ance is denoted as var(7). Here 7 and var(7) are defined 
as 

7 = lim (7^) (10) 

with 
and 

var( 7 ) = lim ((7* ) - ( 7J v) 2 ), (12) 

here Tjv is the transmission rate for system with 2N lay- 
ers (N periods in the corresponding periodic system), and 
notation (• • •) represents the ensemble average. The sam- 
ple size is chosen in such a way that it is much larger than 
the localization length and the ensemble average is car- 
ried out over 200 random configurations. As expected, 
when the randomness is small the LE can be approxi- 
mated by the inverse penertration depth for wave prop- 
agating in the underlying semi-infinite periodic system. 
Like in the case discussed in Ref. Eq], the double peaks 
of var(7) first appear near the gap edges. The peaks of 
var(7) imply the fluctuations of transmission. Further 
increasing A, the var(7) peaks become fatter and flat- 
ter. We observe that if in the vicinity of a peak there is 
another peak, then increasing the randomness will cause 



them to merge. For model 1, as A is increased, a pair 
of var(7) peaks merge with each other inside the gap, 
following the scenario of Ref ]l5p . However, for model 
2 when the randomness is increased the pairs of var(7) 
peaks tend to merge with each other in the pass bands 
and finally destroy the pass bands. It seems that in model 
2 the merging of the double peaks of var(7) in a pass band 
will never be completed. The merging tendency of a pair 
of double peaks merely increase the LE and destroy the 
pass band they belong to. 

We also plot the LE versus its variance in Fig. [| There 
we observe that although for model 1 (with small dielec- 
tric contrast between two layers like that studied in [ 51| ) 
SPS seems still a good approximation in large disorder 
limit, deviation from SPS is clearly observed in model 2. 
Similar quantities have also been calculated for model 3 
and ploted in Fig. || and Fig. S. There we also observe 
novel behaviors of 7 and var(7) and find even larger de- 
viation from SPS as reported in fl7|| . 

From these observations, we find: 

• How does the var(7) vary with the randomness 
depends on many parmeters. For example, the 
impedance constrast, the thickness ratio and which 
kind of layers (for example, high phase velocity, low 
phase velocity etc.) is randomized. 

• When impedance constrast is large, it seems that 
the deviation from SPS is a usual feature. The 
results reported in Jl5[ is based more or less on 
the fact that the impedance contrast between two 
layers is closed to 1. 

In order to understand why the large deviation from 
SPS in model 3 is established we study its transmission 
properties in more detail. In Fig. ^(a) we plot the trans- 
mission rate of the AC waves through a pair of air blocks 
of thickness a. The two air blocks are seperated by a wa- 
ter layer of thickness b — 9999a. Similarly in Fig. ||(b) we 
plot the transmission for the system with 100 air blocks. 
Comparing these two diagrams we find that the air blocks 
are very strong scatterers and thus as few as only two air 
blocks are enough to determine the ranges of band gaps 
and pass bands. If the localization effect of AC waves in 
model 3 is mainly determined by the multiple sacttering 
of AC waves between pairs of air blocks, then one would 
expect that for a system with 100 air blocks the trans- 
mission T/v = ioo in the band gaps can be approximated 
by T]y =2 , here T/v = 2 refers to the transmission rate for 
two air blocks illustrated in ||(a). We indeed observed 
this result in Fig. |^(c). When randomness is not very 
large, this local effect explains why the phase averaging 
process in model 3 is so inefficient that the randomness 
cannot modify the LE much in the band gaps. 
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IV. CONCLUDING REMARKS 

In this paper we studied the statistics of localization 
properties in one-dimensional layered systems. Luapunov 
exponent and its variance are compared for three cho- 
sen models. We find that the band structures of the 
corresponding periodic systems influence the localiza- 
tion properties much if the impedance contrast between 
neighboring layers is not close to 1. In general the single 
parameter scaling is not very accurate and more model- 
dependent parametrs should be included in the detailed 
discriptions of localization behaviors. 

ACKNOWLEDGEMENT 




The work received support from National Science 
Council (No. NSC89-2611-M008-002 and NSC89-2112- 
M008-008). 




0.2 0.4 0.6 0.8 1 1 1.2 1.4 1.6 1.8 2 



Kd/7i k<b>/7t 



P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

V. Baluni and Willemsen, Phys. Rev. A 31, 3358 (1985). 

C. M. Soukoulis, E. N. Economou, G. S. Grest and M. 
H. Cohen, Phys. Rev. Lett. 62, 575 (1989). 

D. Sornette and O. Legrand, J. Acoust. Soc. Am. 92 
(1992). 

A. R. McGurn, K. T. Christensen, F. M. Mueller, and A. 
A. Maradudin, Phys. Rev. B 47, 13120 (1993). 
Z. Ye and A. Alvarez, Phys. Rev. Lett. 80, 3503 (1998). 
Z. Ye, H. Hsu, E. Hoskinson, and A. Alvarez, Chin. J. 
Phys. 37, 343 (1999). 

M. Asch, W. Kohler, G. Papanicolaou, M. Postel, and B. 
White, SIAM Review 33, 519-625 (1991). 
C. H. Hodges and J. Woodhouse, J. Acoust. Soc. Am. 
74, 894 (1983). 

R. Dalichaouch, J. P. Armstrong, S. Schultz, P. M. Platz- 

man and S. L. McCall, Nature, 354, 53 (1991). 

P. W. Anderson, D. J. Thouless, E. Abrahams, and D. 

S. Fisher, Phys. Rev. B 22, 3519 (1980). 

J. Sak and B. Kramer, Phys. Rev. B 24, 1761 (1981). 

J. Flores, P. A. Mello, and G. Monsivais, Phys. Rev. B 

35, 2144 (1987). 

Avraham Cohen, Yehuda Roth, and Boris Shapiro, Phys. 
Rev. B 38, 12125 (1988). 

L. I. Deych, D. Zaslavsky and A. A. Lisyansky, Phys. 
Rev. Lett. 81, 5390 (1998) 

L. I. Deych, A. A. Lisyansky, and B. L. Altshuler, Phys. 
Rev. Lett. 84, 2678 (2000). 

Pi-Gang Luan and Zhen Ye, Phys. Rev. E 63, 066611 
(2001). 



FIG. 1. Band structures (solid lines) and transmission 
curves for the three models disscussed in this paper, (al) 
and (a2) are for the first model, (bl), (b2) are for the second 
model, (cl) and (c2) are for the third model. The broken lines 
in (al)-(cl) represent the inverse penetration depth Im(A'). 
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FIG. 3. var(7) versus 7 for the first and the second models. 



FIG. 2. LE (broken lines) and its variance (solid lines) for 
the first ((A1)-(A5)) and the second ((B1)-(B2)) models. 
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FIG. 4. Lyapunov exponent (broken lines) and its variance 
(solid lines) for the third model. 
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FIG. 6. Transmission vs dimensionless frequency for the 
third model, (a) Transmission for the N=2 system (two air 
blocks), (b) Transmission for the N=100 system (100 air 
blocks), (c) Comparison between T^ =2 (broken line) and 
Tat=ioo (solid line). 
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FIG. 5. var(7) versus 7 for the third model. 
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